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Abstract 

£T) • The inverse scattering transform for a special case of the 3-wave resonant interaction 

equations with non- vanishing boundary conditions is studied. The Jost solutions and the 
fundamental analytic solutions (FAS) for the associated spectral problem are constructed. 
The inverse scattering problem for the Lax operator is formulated as a Riemann-Hilbert 
F ^ , problem on a Riemannian surface. The spectral properties of the Lax operator are for- 

Ph I mulated. 

o\ : 

hh 1 ' 1 Introduction 

: 

One of the important nonlinear models with numerous applications in physics that appeared 
at the early stages of development of the inverse scattering method (ISM), see [35, 23, 36, 28, 
l — *\ 12, 25], is the 3-wave resonant interaction model described by the equations: 

*> '■ , • . * n 

!— + !,,— + xg 2 g 3 = 0, 

CO 1 . • . * n /-i -i \ 

Sn: 1 ar + IVa & +xft « , = ' (L1) 

^t" ! -%J . . dq 3 

o . ^ + + = o. 

Here = qi(x,t), i = 1,2,3, x is the interaction constant, t>j are the group velocities of 
.£h ' the model and the asterisk stays for complex conjugation. The 3-wave equations can be solved 
through the ISM due to the fact that Eq. (1.1) allows a Lax representation (see eq. (2.1) below). 
The main result of the pioneer papers [35] consist also in proving that if Uk(x,t), k = 1,2,3 
satisfy the system (1.1), then the one-parameter family of ordinary differential operators (2.1) 
is iso-spectral (assuming vanishing boundary conditions as \x\ — > oo: ]im.\ x \^ 00 q k (x,t) = 0, 
k = 1,2,3). 

The 3- and iV-wave interaction models describe a special class of wave-wave interactions 
that are not sensitive on the physical nature of the waves and bear an universal character. 
This explains why they find numerous applications in physics and attract the attention of the 
scientific community over the last few decades [35, 36, 9, 23, 24, 25, 7, 8, 10, 29, 11, 18, 4, 5]. 

The interpretation of the inverse scattering method (ISM) as a generalized Fourier trans- 
form and the expansions over the so called squared solutions started in [2] for the nonlinear 
Schrodinger type equations, was soon generalized also for the iV-wave equations [21], see [! j 
and the numerous references therein. It allows one to study all fundamental properties of the 
relevant nonlinear evolutionary equations (NLEEs) which include: 

1. the description of the whole class of NLEE related to a given spectral problem (Lax 
operator L(\) in the form (2.1)) solvable by the ISM; 



1 



2. derivation of the infinite family of integrals of motion associated with L(\); 



3. the Hamiltonian properties of the NLEEs. 

For the case of (1.1) one can show that the model equations are Hamiltonian [35, 36] and 
possess a hierarchy of pairwise compatible Hamiltonian structures [21, 22]. The (canonical) 
Hamiltonian of (1.1) is given by [36, 15, 23, 25]: 

#3-w = \ J°° dx {q k ^ - q* k ^j + K{q,qlql + q^q^ . (1.2) 

Along with the case of (1.1) with vanishing boundary conditions, a special interest deserves 
the case when some (or all) of the functions q k {x,t) tend to a constant as x — > ±oo. Below we 
choose: 

Qi,2(x, t) — > 0, q${x, t) — > pe 19i± , x — > ±oo. 

Here the constants 9 = <p + — </>_ and p are of a physical origin [12] and play a basic role in 
determining the properties of (1.1) with constant boundary conditions and its soliton solutions. 
More specifically, p characterizes the end-points of the continuous spectrum of (1.1) of the Lax 
operator L(\). The discrete spectrum, in this case, may consist of real simple eigenvalues 
k — 1, . . . , N lying in the lacuna — 2p < A/. < 2 p. To them, there correspond the so-called "dark 
solitons" whose properties and behavior substantially differ from the ones of the bright solitons. 
The dark solitons for the nonlinear Schrodinger type equations and their generalizations with 
non- vanishing boundary conditions are studied in [19, 20, 27, 26, 30, 31]. Similar results for 
the discrete nonlinear Schrodinger type equations (the Ablowitz-Ladik hierarchy) are obtained 
in [1,33]. 

It is normal to expect that the properties 1-3, known for the case of vanishing boundary 
conditions will have their counterparts for the case of constant boundary conditions. However 
there is no easy and direct way to do so. The most concise and systematic treatment of both 
problems (on the example of nonlinear Schrodinger equation) is given in [12]. It is shown there 
that one may relate both cases by taking a limit p — > 0. Of course, in this limit most of the 
difficulties, related mostly with the end-points of the continuous spectrum disappear. From 
[12] one can see that the spectral data, the analyticity properties of the Jost solutions and 
the corresponding Riemann-Hilbert [13, 34, 32] problem are substantially different and more 
difficult for p > 0. 

The aim of the present paper is to study the direct scattering problem for the Lax operator 
and its spectral properties. In Section 2 we start with the Lax representation and the construc- 
tion of the Jost solutions of the Lax operator L. In Section 3 we outline the construction of the 
fundamental analytic solutions (FAS) of L. We also formulate the Riemann-Hilbert problem 
on the relevant Riemannian surface, that is satisfied by the FAS. In section 4 we derive the 
time evolution for the scattering matrix. Section 5 is devoted to constructing the resolvent of 
L in terms of the FAS and the spectral properties of L. The effects of the boundary conditions 
on the conserved quantities of the 3- wave equations are analyzed in Section 6. We finish with 
a brief discussion and conclusions. 



2 Lax representation and Jost solutions 

The idea of the Inverse Scattering Method (ISM) is based on the possibility to linearise the 
nonlinear evolutionary equation (NLEE) [3, 6, 36]. To this end we consider the solution of the 
NLEE Uk(x, t), k = 1, 2, 3 as a potential in the Lax operator L(\). 
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Consider the pair of Lax operators: 

14 = {£- + [J, Q{x, t)} - XJj ^{x, t, A) = 0, (2.1) 
Mip ee + [/, Q(x, t)} - Xlj ij(x, t, A) = 0, 

with 

n (°* q n q3 \ J = diag(J 1 ,J 2 ,J 3 ), 

Q = ?i 92 , r , r T T v (2.2) 

/ / = diag(/ 1 ,J 2 ,/ 3 ). 

We presume that the potential matrices Q(x,t), I and J are traceless (i.e. the Lax operators 
take values in the algebra sl(3, C)) and the eigenvalues of J and J are ordered as follows: 
Ji > J 2 > J3, h > I 2 > h (Ji + J2 + J3 = and ^ + J 2 + 1 3 = 0). Here A G C is a spectral 
parameter. 

The compatibility condition for (2.1) leads to 

i[J, Q t ] - Q x ] + [[/, Q], [J, Q}] = 0, (2.3) 

which is equivalent to (1.1), if the potential matrices are taken from (2.2). The group velocities 
in (1.1) take the form: 

h — h h — h h — h 



Vl = -? D v 2 = ~? T~i v 3 



J\ — J2 J2 — J3 J\ — J3 

while the interaction constant x reads: 

x = JJ 2 + J 2 h + J3I1 ~ -hh - -hh - Jih- 

We will assume also, that the potential of the Lax operator is a subject of constant boundary 
conditions as \x\ — >■ 00: 

lim qi{x,t) = lim q 2 (x,t) = 0, lim q 3 (x,t) = qf = pe 1 ^. (2.4) 

x— >±oo x— >±oo x~>±oo 

Equivalently, for the potential matrix Q(x,t) one can write: 



lim Q(x,t) =Q±, Q±=\ 

x— >±oo 



pe^± 



pe -i<f>± 



The difference 9 = <p + — 0_ of the asymptotic phases <fi± plays a crucial role in the Hamiltonian 
formulation of the 3-wave model with constant boundary conditions: its values label the leaf 
on the phase space M of the model (1.1) where one can determine the class of admissible 
functionals, and to construct a Hamiltonian formulation. The two asymptotic potentials Q± 
are related by 

Q+ = Q{e)Q-{t)Q-\e), (2.5) 

where 9 = <p + — 0_ and 

■ e ie/2 Q 

Q(9) =| 1 I . (2.6) 
0e- w / 2 



The direct and the inverse scattering problem for the Lax operator (2.4) will be done for fixed 
t and in most of the corresponding formulae t will be omitted. 

The starting point in developing the direct scattering transform for the Lax operator (2.1) 
are the eigenf unctions (the so-called Jost solutions) of the auxiliary spectral problem 



L(x, t, \)ip±(x, t, A) = 0, 
determined uniquely by their asymptotic behavior for x — > ±oo respectively: 

lim ^±M,A)e iJ M* = V; ± ,o(A)P(A), 

x— >±oo 

where -P(A) is a projector: 

P(A) = diag(0(|ReA| -2p),l,0(|ReA| -2p)) 



(2.7) 



{2.1 



(2.9) 



and 9(z) is the step function. As we shall see below, -P(A) ensures that the continuous spectrum 
of L has multiplicity 3 for |Re A| — 2p > and multiplicity 1, for — 2p < Re A < 2p. The x and 
t-independent matrices ip±fi{ty m (2-8) diagonalize the asymptotic Lax operators: 



d 

L±{x,t,\) = i— + [J,Q±] 



A J. 



Indeed, 



where 



{[J,Q±\ - AJ)^±,o(A) = -V>±,o(A)J(A), 



(2.10) 



(2.11) 



J(A) = -diag(J 1 (A),J 2 (A),J 3 (A)), 

JiW = \ [^A + (Ji - J 3 W\*-4p? 
1 



J2(A) 



-AJ 2 



(2.12) 



J 3 (A) = - [J 2 A - (J x - J 3 ) v / A 2 -4p 2 



jfe(A) = VA 2 - 4p 2 



For the choice of Q± as in (2.5) we have: 

1 



2p 0-(A + fc)e^± 
^±o(A) = ^ x I 1 

V 2X (k + A) \ (A + jfe) e -**± 2p 



(2.13) 



Here and below we will deal with the Riemannian surface related to /c(A); its first sheet is fixed 
up by the condition: signlmfc(A) = signlmA. The Jost solutions tp + (x,t, A) and ip-(x,t, A) are 
related by the scattering matrix T(t, A): 



T(t, A) = ip~, l {x,t, \)ip_(x,t, A), detT(A) = 1. 



(2.14) 



For a sake of convenience, from now on, instead of the spectral parameter A we will be using 
the so-called " uniformizing variable" 



C = -(A + A;(A)). 



(2.15) 



In terms of ( we have: 



A = P (C + ^ 



fc(A)=p(C-j 
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Then, the formulas (2.12) take the form: 

J(C) = -pdiag (.M + j, J2 (^)' JlC+ f)' (2 - 16) 
Along with the Jost solutions (2.8) it is convenient to consider slightly modified Jost solutions: 
V±(xX) = ^o(OM^C)e i - mx , lim »7±foC) = l, (2.17) 

and satisfying the following associated to (2.7) equation: 

i^- + ViifJ, Q(x, t) - Q ± }^ ±fi V± ~ P(C + C 1 )[J(0,V±(x, t, ()] = 0. (2.18) 

Equivalently, f]±{x, () can be regarded as solutions to the following Volterra-type integral equa- 
tions: 



v±(x,C) = a + i / dyc- iJ 

J ±00 

x 1>£ [J, Q(y) - Q±]V±,o77±(2/, ()e iJ ^ x - y \ ( e R, (2.19) 



where the diagonal matrix J(£) is given in (2.16). 

In addition the second column of the Jost solution is defined also on the unit circle in the 
C-plane: 



,(2)v /~\ — x, i_ /)' r ^,„-i(^(C)-J2(C))(x-!/) 



(r)±)k2(x,C) = h2 + i dye' 

J±oo 

x teo[ J > Q(y) - Q±]^±,o^±(j/, 0) , ICI = l. (2.20) 
3 The fundamental analytic solutions of L 

In order to construct the fundamental analytic solutions (FAS) of L we first need to determine 
the regions of the complex (-plane in which the imaginary parts of the eigenvalues of J(() are 
ordered. To do this we first need to find the curves on which 

Im(J,(C) - MO) = 0, 1 < J < A; < 3; ImJ^C) = 0. (3.1) 

Writing down ( = Kle*^ we find: 

t t \ i ai , d\ — J2 



MJi(C) 


- MO) 


Im(J 2 (C) 


- MO) 


Im(Ji(C) 


- MO) 




ImJ 2 (C) 



ICI 

J2-J3 

ICI 

1 



sm <p , 
sin 0o, 



(3.2) 



'"ICI 

-J2P (jCl - j^J sin 0o- 

Since J\ — J 2 > 0, J\ — J3 > and J 2 — J 3 > it is easy to see that the solutions of the eqs. 
(3.1) are O — and O = ^ i- e - the rea l ax is m the complex (-plane; in addition ImJ 2 (() = 
for ((I 2 = 1. 
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The complex (-plane is split into four regions f^, k — 1, . . . , 4 formed by the intersections 
of the upper and lower complex half-planes C + and C_ with the unit circle S, see figure 1. The 
ordering of ImJfc(A) in each of them depends on the sign of J 2 and are as follows: 



^4 



ImJi(A) > > ImJ 2 (A) > ImJ 3 (A), 
ImJi(A) > ImJ 2 (A) > > ImJ 3 (A), 
ImJ 3 (A) > > ImJ 2 (A) > ImJ^A), 
ImJ 3 (A) > ImJ 2 (A) > > ImJi(A), 



(3.3) 



for J 2 > and 



(3.4) 



fii : ImJi(A) > ImJ 2 (A) > > ImJ 3 (A), 

tt 2 : ImJi(A) > > ImJ 2 (A) > ImJ 3 (A), 

fi 3 : ImJ 3 (A) > W 2 (A) > > ImJ^A), 

fi 4 : ImJ 3 (A) > > ImJ 2 (A) > ImJi(A), 

for J 2 < 0. 

Let us first construct the FAS in the region f^. Following the ideas of [36, 22] we introduce 
it as the solution of the following set of integral equations: 



6+(x,C)j =5 kl + i f dye -WyJi(0)(*-y) 

)kl Joe 

x U^[J,Q(y)-Q + ]^ +t o^ 1) (y,0 



kl 



k<l. (3.5) 



-*(J*K)-Ji(0)(*-w) 



k>i. 



kl 



The proof of the fact that £^{x, C) is an analytic function of ( for any £ G Qi is based on the 
fact, that due the ordering (3.3) all exponential factors in eqs. (3.5) for ( G C + are decaying. 
This ensures the convergence of all integrals as well as the existence of Cm( x ^ f° r a ^ C e 
Similar conclusions can be drawn also for the (d/d() k £,ty Indeed, taking derivatives of (3.5) 
with respect to ( would give rise to terms polynomial in x in the integrands. Such terms, 
however, are suppressed by the exponential factors, which allows to conclude, that along with 
£m( x iC) a l so its derivatives (d/d() k £,^ exist for any positive k. Thus, we briefly outlined the 
idea of the proof that Xn){ x i C) — C) e are FAS of the Lax operator L for £ G fii- 

Similarly, the FAS in the region U4 is the solution of the set of integral equations: 



£(4)( x '0 L = Sia + i I dye 



ki 



-t(Jfc(0-Jj(C))(*-w) 



x \^[J,Q(y)-Q-]^ + , ^ 4) (yX) 



kl 



k<l. (3.6) 



{f(4)foC) 



kl 



dye 



-i(J h (C)-Ji{C))(x-y) 



k > I. 



ki 
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The proof of the analyticity of £u\(x, C) f° r an Y C £ ^4 is similar to the one for £X\(x, ()■ 

It remains to outline the construction of £m(x, A) and £Xs(x, A) for the regions S7 2 and Jl 3 . 
To this end we make use of the involution of the Lax operator L that is a consequence of 
Q = —Q>. Then we conclude: 

Xp)(s,0 = (x#rW/O, X$)(*,C) =(x (4 ' ) V 1 (x,l/C). (3.7) 

The next step is to analyze the interrelations between the Jost solutions ip±(x, and the FAS 
X + {x, C) an d X - !^) 0- It is natural to expect that they are linearly related. Skipping the details 
we note that: 

= <M^0%(0, 4^,0 = M*, tK a) (()D+ a) ((), 

X w (x,() =^-(x,()s m ((), xlp ) (x,0 = Mx,()T$ )) (()D(p ) ((), 

where (gRU§, a=l,3 and — 2, 4 match the indices of the regions of analyticity flk- We 
will often omit the indices a and /3, since their values are clear from the figure 1. Here and 
(resp. STg) and TZA are upper triangular (resp. lower triangular) matrices whose diagonal 
elements are all equal to 1; the matrices D,s, D,^ are diagonal ones. In fact these matrices 
are directly related to the Gauss decomposition of the scattering matrix T(0 (2.14): 



T{C)=T*{QD ± {Q{S ± {())-\ (3.9) 



where: 



1T 1 + (C)T 3 + (C)\ / 1 0' 

T + (C) = | 1 T+(C) , T~(t,0 = Tf(C) 1 

I J V T 3~(C)T 2 (C)1, 

ist(C)si(0\ / 1 oo' 

5 + (o= I o i s+(o , ^-(t,c)= sr(o i o 

oo i J \5 S -(C) fliCC) l 

™ 2 + (0 1 



(3.10) 



D + (C) = diag m+(C), 



(3.11) 



Here mf(() are the principal upper/lower minors of order k of the scattering matrix (2.14); the 
explicit expressions for the matrix elements of and S ± in terms of Tij(() are given by: 

q+(n = t 12 (Q = ^(0^(0-^(0^(0 

llU T n (0' 21 J ^(0^(0-^(0^(0' 

T 12 (QT 23 (0-r 13 (OT 22 (0 

3lu r n (0T 22 (C)-r 12 (or 2 i(C)' 1 ' ' 

{ ) = 721(0 T - (c) _ ^(OTuCC) - r 3 i(C)Ti 2 (0 



Tn(0' 2 ' 7 T n (0T 22 (0 -T 12 (OT 21 (0' 



T s(0 = ^77T. ( 3 - 13 ) 
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and 



o- rn _ r 33(C)r 21 (C)-T 2 3(C)r 31 (C) 
llu T 32 (C)T 2 3(C)-T33(C)r 22 (C)' 1 • ) 



q-tn = TM <r(n = t 3i(O t ^(0-t 32 (C)t 21 (c) 

2lU " T 33 (C)' 3lU T3 2 (C)T 23 (C)-T33(C)T 22 (C)' 

T - (n = T 12 (C)T33(C)-T 13 (C)T3 2 (C) , , 

1 lu rii(C)r 3 3(C)-T 13 (C)r3i(C)' 1 ' ] 

2{Q ~T 33 (cr 3(g "T 33 (o' 

Due to the the special choice of the matrix Q(x,t) (2.2) it follows that S~(Q = (S + (1/(*)Y 
and T + (() = (T~(l/£*))t, so the matrix elements of S~(() and T + (() can be reconstructed 
from (3.12). 

One of the most effective method for solving the inverse scattering problem for a given Lax 
operator L is to reduce it to a RHP [ ]. On the complex (-plane it can be formulated as 
follows: 

G(() = e -iAO^-inOtfs~)~ 1 S + e iJ ^ x+iF ^ t . 

The relation (3.16) holds true for k £ R in the complex fc-plane. The RHP for Lax operators 
with vanishing boundary conditions look similarly. However, the relation (3.16) is more compli- 
cated due to the fact, that we are dealing with an RHP formulated on the Riemannian surface 
related to the root k(X) = a/A 2 — 4p 2 . 

The sewing function G(x, () gives the minimal set of scattering data, sufficient to reconstruct 
the scattering matrix T(£). 

One of the important uses of the RHP is that it allows one to use the Zakharov-Shabat 
dressing method and construct the soliton solutions of the relevant NLEE. Doing this one 
should specify the dressing factor as a rational function of the uniformizing variable ( rather 
than A. 



4 The time evolution of the scattering matrix 

We start by noting that we could use a bit more general M-operator than the one in (2.1), 
namely: 

Mip = (i^ + [I, Q(x, t)} - C/) t, = 1>{x, t, C)F(C). (4.1) 

The compatibility condition [L, M] = holds true for any x- and t-independent matrix F{Q. 
We will fix up F(Q, requiring that the definition of the Jost solutions (2.8) holds true for all 
t. Let us identify ip(x,t,Q as ip + (x,t,() (resp. ip-(x,t,Q) and take the limit x — > oo (resp. 
x — > — oo). This gives: 

(iI,Q±}-W±,o = Tp±,oF((). (4.2) 

It is easy to check that ip± y o diagonalize also [I, Q±] — (I and therefore F(() is a diagonal 
matrix: 

F(C) = diag(/i(C),/ a (C),/s(0), 
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where fj(£) are the eigenvalues of [/, Q±] — (I. In terms of ( we have: 



F(C) = pdiag [h(+^, J 3C + ^)- (4-3) 



Next, we insert ^ = ift+(x,tX) into (4.1) and take the limit x — > — oo. Thus we obtain, that 
if the matrix elements Uk{x,t), k = 1,2, 3 of the potential of L(() satisfy the NLEE (1.1) then 
the time evolution of the associated scattering matrix is given by the linear ODE: 

i^-[n<),T{t,O] = 0- (4-4) 
As a consequence the Gauss factors of T(t, () satisfy 

<^-[F(c),r t (t,C)] = o, i^-[F(Q i S±(t i Q] = 0, ^ =0, (4.5) 

From the last equation it follows that the principle minors m^(Q and mf(Q of the scattering 
matrix are time-independent and can be considered as generating functionals of the integrals 
of motion for (1.1), while for the off-diagonal ones we get: 

T^t, = 1^(0, C) e -*(*(0-/i(0)f . (4. 6 ) 

The function F(() is known as the dispersion law for the 3- wave equations with constant 
boundary conditions. 



5 Spectral Properties of the Lax Operator 

The crucial fact that determines the spectral properties of the operator L(() is the choice of the 
class of functions where from we shall choose the potential Q(x). Below, for a sake of simplicity, 
we assume that Q(x,t) satisfies 

Condition C.l Q(x,t) is smooth for all x and t and is such that 

lim \x\ p (Q(x,t) -Q±) = for all p = 0,l,.... 

x— >±oo 



The FAS x ± ( a; ! C) °f -^(C) allows one to construct the resolvent of the operator L and then 
to investigate its spectral properties. By a resolvent of L(() we understand an integral operator 
R(C) with kernel R(x, y, () which satisfies 

L(0(R(C)f)(x) = f(x), (5.1) 

where f(x) is an 3-component vector complex- valued function with bounded norm, i.e. 

/oo 
dy\f T (y)f(y)\ <oo. 
-oo 

From the general theory of linear operators we know that the point ( in the complex £-plane 
is a regular point if -R(C) is a bounded integral operator. In each connected subset of regular 
points -R(C) is analytic in (. The points ( which are not regular constitute the spectrum of 
L(C). Roughly speaking, the spectrum of L(() consist of two types of points: 
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© 



© 



c. 




-2p 



2p 



Figure 1: The continuous spectrum of L in the complex A plane (left panel) and in the complex 
C plane (right panel), 

• i) the continuous spectrum of L(() consists of all points ( for which R(() is an unbounded 
integral operator; 

• ii) the discrete spectrum of L(() consists of all points ( for which R(() develops pole 
singularities. 

Let us now show how the resolvent R(() can be expressed through the FAS of L((). Indeed, 
if we write down R(() in the form: 




(5.2) 



the kernel R(x, y, () of the resolvent is given by: 



R(a) (X,V,C) = R (a) ( x ,y,0, for c G fi (a) • 



(5.3) 



where 



t)(^2/.C) = -ixf a) (a:,C)©f a )(^-S/)(xf Q )) "(vX), C^(a), 



(5.4) 



and 



Q+(x -y) = diag (-% - x), 6(x - y), 9{x - y)), 
= di&g(-6(y-x),-6(y-x),6(x-y)), 

®t( x -y) = dia s (°( x - v), -°(y - x ), -°(y - ^))> 

&2(x-y) = diag (-9(x - y),9(x - y), -9(y - x)), 



(5.5) 



for J 2 > 0, and 



Q+(x -y) = diag {-9{y - x), -9(y - x), 9{x - y)), 

©2 ( x -y) = dia s (-% - x ), °( x - y), °( x - v)), 

®t( x ~y) = diag (9(x - y), 9(x - y), -9(y - x)), 



(5.6) 



@ 4 (x-y) = diag {-9(x - y), -9(y - x), -9{y - x)), 
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for J2 < 0. The next theorem establishes that R(x, y, () is indeed the kernel of the resolvent of 

Theorem 1. Let Q(x) satisfy the conditions (C.l) and is such that the minors ^i^(C) have 
a finite number of simple zeroes (p. Then 

1. R± (x,y,C) is an analytic function of ( for (gC ± having pole singularities at £■ ; 

2. R (x, y, is a kernel of a bounded integral operator for ( G R U £; 

3. R(x, y, () is uniformly bounded function for ( G RU £ and provides a kernel of an un- 
bounded integral operator; 

4- R ± (x,y,() satisfy the equation: 

L(C)R ± (x,y,C) = *5{x-y). (5.7) 

Proof: 

1. is obvious from the fact that x ± ( x , C) are the FAS of L((); 

2. Assume that ( G fii and consider the asymptotic behavior of R + (x,y,() for x,y — >■ 00. 
From equations (3.8) and (5.4) we find that 

3 

p=l 

Due to the fact that x + (x, () has triangular asymptotics for x — > 00 and ( G C + and for 
the correct choice of <d + (x — y) (5.5) we check that the right hand side of (5.8) falls off 
exponentially for x — > 00 and an arbitrary choice of y. All other possibilities are treated 
analogously [14, 21]. 

3. For ( G tU§ the arguments of 2) can not be applied because the exponentials in the 
right hand side of (5.8) only oscillate. Thus we conclude that R ± {x ) y ) Q for £ G K U S 
is only a bounded function and thus the corresponding operator -R(C) is an unbounded 
integral operator. 

4. The proof of eq. (5.7) follows from the fact that L(Qx ± ( x ^ C) = an d 

= «(*-»). (5.8) 

ax 



Thus we conclude that the continuous spectrum of L in the complex /c-plane coincides with 
the contour of the RHP MUS with multiplicity 3 on R and multiplicity 1 on S. On the complex 
A-plane the continuous spectrum of L is on the real axis; it has multiplicity 3 on the semi-axis 
Re A < — 2p and Re A > 2p and multiplicity 1 in the 'lacuna' — 2p < Re A < 2p. 
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6 Conserved Quantities for the 3-wave Equations 



As we already mentioned above, the diagonal factors D ± (() are time independent and can 
be used to generate the infinite set of integrals of motion for (1.1). For the 3-wave resonant 
interaction equations, these matrices are expressed through the principal upper/lower minors 
m ± (C) of the scattering matrix T(£) (2.14). Skipping the details (see [21]) we get: 



]nD± A = -~{J k - J k+1 )? k + {Jt - J 3 )? 3 , (6.1) 
The momenta k — 1, 2, 3 are given by: 

/OO /*00 POO 

dx\ qi (x)\\ T 2 = dx\q 2 (x)\ 2 , 0> 3 = / dx{\q 3 {x)\ 2 -p 2 ). (6.2) 
-oo J— oo J— oo 

The fact that In mf generates integrals of motion can be considered as natural analog of the 
Manley-Rowe relations [35, 23]. In the case of (2.2), then (6.1) is equivalent to the existence 
of two additional first integrals for the model (1.1) 

h = {Ji ~ J2P1 + (Ji ~ J3P3 = const, 

h = (J2 - J3P2 + (Ji - Js)7 3 = const, 1 ' ' 

which are linear combinations of the momenta (6.2) and can be interpreted as relations between 
the densities \q a \ 2 of the waves of type a. The total momentum for the 3-waves is also a 
conserved quantity: 

V={Jx- J2P1 + (J 2 - ^3)^2 + (Ji - -h)% = const. (6.4) 

The integral of motion D 2 is proportional to the Hamiltonian of the 3-wave equations (1.2). 

For the case of constant boundary conditions, the functional CP3 is subject of regularization 
(6.1) by using the asymptotic values of the potential Q±, while the functional H 3 _ w remains 
the same as for the case of vanishing boundary conditions. 



7 Conclusions 

We studied the direct scattering problem for the Lax operator and its spectral properties. This 
includes: the construction of Lax representation and the Jost solutions of the Lax operator L. 
Furthermore, we outlined the construction of the fundamental analytic solutions (FAS) of L 
and formulated a Riemann-Hilbert problem for the FAS on a relevant Riemannian surface. 

We also outlined the construction of the resolvent of L(() in terms of the FAS and the 
spectral properties of L. Finally, we briefly discuss the effects of the boundary conditions on 
the conserved quantities of the 3-wave equations: we showed that the total momentum for 
non-vanishing boundary conditions needs regularization, while the Hamiltonian remains the 
same. 

Similar analysis can be done for a 3-wave resonant interaction model with more general 
boundary conditions: lim^^.oo q k (x, t) = q k (k = 1 or 2) and \im x ^ +CJO q 3 (x,t) = q 3 . This may 
require the matrices Q(9) (2.6) to have also off-diagonal entries. 

It is an open problem to to derive the soliton solutions of (1.1) in the case of constant 
boundary conditions (the so-called "dark solitons"), the dark-dark and dark-bright soliton 
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solutions [30] by modifying the dressing Zakharov-Shabat method [ ], or by using the Darboux 
transformation method [10]. 

Another challenge is to extend this analysis also for systems, describing resonant interactions 
of iV waves [36, 16] or to A-wave type systems related to simple Lie algebras [16, 17]. 

Another open problem is to study the behavior of the scattering data at the end-points of 
the continuous spectrum in the complex A-plane; this requires generalization of the the method 
developed in [12]. 
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